Subleading shape functions in B — > X u £u and the determination of \V u b\ 
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We calculate the subleading twist contributions to the endpoint of the inclusive lepton energy 
spectrum in B — > X u tv. We show that the same two subleading twist functions that appear in 
the decay B — > X s y govern the subleading effects in this decay. Using these results we find large 
0(Aqcd/?tis) corrections to the determination of \V u b\ from the endpoint of the charged lepton 
spectrum. Using a simple model for the relevant subleading shape functions, we estimate the 
uncertainty in \V u b\ from Aqcd /m& corrections to be at the ~ 15% level for a lower lepton energy 
cut of 2.2 GeV. 



A precise determination of the CKM matrix element 
\V ub \ is one of the most important measurements in 
the current quest to overconstrain the unitarity tri- 
angle to test the standard model. The theoretically 
cleanest determination of \V u b\ comes from the inclu- 
sive rate for B — > X u lv, which may be calculated 
model-independently via an operator product expansion 
in terms of local operators Such a measurement is 
hampered by the experimental necessity to suppress the 
overwhelming charm background, which usually requires 
tight kinematic cuts. With a few exceptions Q], such 
cuts invalidate the local OPE usually used to calculate 
the inclusive branching ratio, and a twist expansion has 
to be performed instead ^ . At leading order in this twist 
expansion, and ignoring perturbative corrections, the dif- 
ferential lepton energy spectrum in inclusive B — > X u £u 
decay is given by 
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where f(w) is the nonperturbative light cone structure 
function of the B meson 
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Since this function is a property of the B meson, it 
is process independent and can thus be measured else- 
where. The best way to measure this structure function 
is from the photon energy spectrum of the inclusive de- 
cay B — > X s j 0. Up to perturbative and subleading 
twist corrections, this spectrum is directly proportional 
to the structure function, 
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Thus, combining data on B —> X s j with data from 
B — > X u lv, one can eliminate the dependence on the 
structure function and therefore determine \y ub \ with no 
model dependence at leading order j| |ij |f| . The relation 
is 
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Recently, the CLEO collaboration M measured the rate 
for B — » X u lv decay in the endpoint interval 2.2 GeV < 
Eg < 2.6 GeV and combined this with their measurement 
of the B — > X s j photon spectrum to determine the value 

|Kb| = (4.08 ± 0.34 ± 0.44 ± 0.16 ± 0.24) x 10" 3 (6) 

where the first two uncertainties are experimental, and 
the third corresponds to the theoretical uncertainty in 
the relation between T(B — > X u lv) and \y ub \. The fourth 
uncertainty is an estimate of the nonperturbative uncer- 
tainties in the relation (Q), and is determined by vary- 
ing the parameters of f(u>) by 10% (as expected for an 
0(AQCD/wfc) correction). 

In this paper we calculate the subleading twist cor- 
rections to Eq. (|4). While perturbative corrections to 
this relation have been studied in detail ||, || , much less 
is known about nonperturbative corrections suppressed 
by powers of Aqcd/?™&- In @ we introduced a formal- 
ism (similar to that used for higher twist calculations 
in deep inelastic scattering |^|) which allows these sub- 
leading twist corrections to be parameterized in terms 
of subleading shape functions, and calculated the sub- 
leading twist corrections to the photon energy spectrum 
in B — > A, 7. These corrections are suppressed by one 
power of Aqcd /fn b and can be written in terms of three 
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FIG. 1: The two full theory graphs required for the matching 
onto the subleading shape functions 
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and the functions F, h\ and H2 will be defined later. 

The subleading twist corrections to the lepton energy 
spectrum in B — + X u lv may be calculated in much 
the same way as the photon spectrum for B — > X s ~/. 
The matching calculation is performed by expanding full 
QCD matrix elements to subleading twist and identify- 
ing the result with non-local operators. The operator 
mediating the decay B — » X u iv is given by 



O = {ub)v-A{£v) V -A = {uv) V -A{ib)v-A , 
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where we have used a Fierz transformation to obtain the 
right hand side of Eq. (g). Using the optical theorem, 
the differential lepton energy spectrum can be obtained 
from the imaginary part of the forward scattering matrix 
element 
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The phase space integrals are automatically performed 
by taking the imaginary part of one loop diagram shown 
in in Fig. |l|(a). Calculating this diagram and taking the 
imaginary part we find 
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and we further defined p^ = — p^. Finally, con- 
tains the contraction of the lepton spinors 

Lu V = j^Pl i> g IvPl , (13) 
with P L = (1 -7 5 )/2. 

To expand this expression to subleading twist, we iden- 
tify the scaling of the various components of the mo- 
menta. Defining 



pi = m b V^ + , Q>* = m b v>* - p^ . 



(14) 



as well as the light cone vectors n M , n M (n 2 = n 2 = 0, 
n ■ n = 2, 2v — n + n) by 



[hi(m b - 2E 1 ) + H 2 {m b - 2JS r )](7) we find 
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fc" ~ 0(A QCD ) 

(this should be compared with the power counting for 
the local OPE, in which n ■ Q scales as m b ). Using this 
momentum scaling in Eq. ( |TT| ) and expanding in powers 
of Aqcd /m b we find at leading order 
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For pseudoscalar B meson decays, only the parity even 
part of this amplitude contributes. The matching onto 
the low energy theory is obtained by comparing this with 
the general expression at leading twist 
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where P + = (1 + f)/2 and 

O (u>) = hS(in ■ D + ui)h. (19) 

Thus, we find 

C a {w) = 2T 6{m h -2E t -u), (20) 

in agreement with Eq. (|l|). 

At subleading twist there are five nonlocal operators 
contributing to B meson decays |7j 
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where s a = P+7 a 7 D P+. The Feynman rules for these operators were given in n ■ A = gauge in Ref. § (note that 
we have used slightly different definitions of O3 and P4 from Ref. 0). To match onto these operators we expand the 
full theory calculations to one order higher in k^/uib and n ■ Q/rrib. For the graph with no external gluon we find 
from (IrT 
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Of the above operators only 0^{w) and 6)3(0;) contribute 
to the zero gluon matrix element. In addition, we use the 
spinor expansion relating the full QCD spinors u to the 
two component HQET spinors h 



u= 1 



2mb 



(24) 



and as before only the parity even part of this amplitude 
contributes to the decay rate. This allows us to write 



Cf(w) = ^-ni*d(mb - 2En - w) (26) 
2F 

C 3 M = -9(m b - 2E e - u) . (27) 



The result for the the coefficient 63(0;) agrees with the 
prediction from reparameterization invariance, which re- 
lates this coefficient to the leading order Wilson coeffi- 
cient to all orders in perturbation theory 
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+k 2 ± 6{n-k + m b -2Ei)] . 

From this we can easily read off the Wilson coefficients 
for the subleading twist operators 0^{uS) and 03(0*). We 



C 3 (u) = -C (u) . (28) 
To determine the matching coefficients of the opera- 
tors P^ a (o;) and P"(ui), we need to calculate the matrix 
element with an explicit gluon in the final state. The 
required diagram is shown in Fig. 0(b). Performing the 
calculation in n ■ A = gauge and expanding to sublead- 
ing twist gives 
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The first and second term in this expression correctly 
reproduce the one gluon terms from the operators 0% (u>) 
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The Wilson coefficients for the operators 
P^ q (cl>) and P"(w) are determined from this expression 
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To calculate the differential decay rate to subleading 
twist, we need expressions for the matrix elements of 
these operators. These matrix elements may be param- 
eterized in terms of four subleading structure functions 
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There are some additional simplifications that may be 
used to reduce the number of unknown functions. First, 
the function t(uj) comes from expanding the Lagrangian 
to higher order in 1/uib and therefore always occurs in 
the same linear combination f(tu) + t(oj) /2mb- Together 
with the reparametrization invariance constraint (^8|) , it 
therefore follows that for any B meson decay one always 
finds the same linear combination of leading and sublead- 
ing shape functions 

F(w) = /H + -L [t(«) - 2 G 2 (w)] . (34) 

Putting all this information together, we therefore find 
for the lepton energy spectrum in B — > X U £D 
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Eq. (|||) is the central result of this paper. We can check 
this result against the known result from the local OPE 
by using the expansion of the nonperturbative functions 
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H 2 (lj) = -X 2 S'{uj) + ... , (36) 
where . . . denote terms 0(Aq CD /m 2 ) in the twist expan- 
sion. This gives 
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Here 3r = —2(71 + 37^) and the pi, % are the papame- 
ters used in jlj, [p| . This agrees with the corresponding 
expressions from Rcfs. O, p"2] ]. 

We can now use these results to get the 0(Aqcd/wb) 
corrections to the determination of \V u b\ from the com- 
bined measurement of the semileptonic and radiative B 
decays given in Eq. (^) . Using the differential decay rates 
given in Eqs. ([?]) and (p5|), together with Eq. (JsJ) , we find 
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where rg and To are defined in Eqs. (||) and (|12j), respec- 
tively. This finally leads to 
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FIG. 2: The subleading twist corrections to the \V u b\ relation 
( |i~lj ) as a function on the lepton energy cut, using the simple 
model described in the text. The vertical line corresponds to 
the kinematical endpoint of the semileptonic b — > c decay. 

where S(E C ) is the correction to the extraction of \V u t\ 
from subleading twist. 

We note from Eq. (|36| ) that the first moments of hi(u>) 
and H2(u>) have opposite signs, and so we expect their 
effects to add, not cancel, in the expression ([4l]). This 
gives the naive estimate S ~ 2Aqcd /mb ~ 0.2, so the 



contributions of the subleading twist shape functions are 
sizable. 

To get a more quantitative estimate for these correc- 
tions, we can use the simple model for the subleading 
structure functions introduced in |Q , in which these func- 
tions are taken as the derivative of the leading order 
shape function, normalized such that they reproduce the 
known first moments. Using the one parameter model 
for the leading order structure function given in |fl4| 
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(where A = ms — mj,) this gives 
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and leads to the expression for the subleading twist con- 
tributions 
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Using the numerical value A = 0.47 GeV for the model 
parameter, we can estimate the subleading twist con- 
tributions 8 as a function of the lepton energy cut E c . 
This is shown in Fig. [|. For the cut used by the CLEO 
collaboration || to reject background from semileptonic 
b — > c transitions E c = 2.2 GeV, the subleading twist 
corrections to \V u b\ are about 15%. Since this number is 
strongly model dependent, we estimate the uncertainties 
in \V u b\ to also be at the ~ 15% level. From Fig. || one 
can see that these uncertainties grow rapidly if the cut 
is raised further. Conversely, lowering the cut below 2.2 
GeV reduces the uncertainty from subleading twist con- 
tributions. However, as the cut is lowered below 2.2 GeV 
the effect of the resummation of local operators in the 
twist expansion becomes small and the usual local OPE 
in terms of local operators becomes valid, eliminating the 
need to extract f(w) from the B — » X s j spectrum. 



We would like to thank the staff of the Aspen Center 
for Physics for their hospitality while some of this work 
was completed. This work was supported by the Depart- 
ment of Energy under grants DOE-FG03-97ER40546 and 
DOE-ER-40682-143 and the Natural Sciences and Engi- 
neering Research Council of Canada, and by the DFG 
Graduiertenkolleg "Hochcnergiephysik and Teilchenas- 
trophysik" , from the DFG Forschergruppe "Quanten- 
feldtheorie, Computeralgebra und Monte Carlo Simu- 
lationen" and from the Ministerium fur Bildung und 
Forschung bmb+f. 

Note: While this work was being completed, we became 
aware of similar results by Leibovich, Ligeti and Wise 
J7^/. We are grateful to these authors for sharing their 
results with us before publication. 



[1] J. Chay, H. Georgi and B. Grinstein, Phys. Lett. B247 M.A. Shifman and M.B. Voloshin, Sov. J. Nucl. Phys. 41 

(1990) 399; 



6 



(1985) 120; 

LI. Bigi et al, Phys. Lett. B293 (1992) 430 [Erratum- [8] 
ibid. B 297, 430 (1992)]; 

LI. Bigi et al, Phys. Rev. Lett. 71 (1993) 496; [9] 
A.V. Manohar and M.B. Wise, Phys. Rev. D49 (1994) [10] 
1310. 

[2] C. W. Bauer, Z. Ligeti and M. E. Luke, Phys. Lett. B [11] 
479, 395 (2000); Phys. Rev. D 64, 113004 (2001) 

[3] I. I. Bigi, M. A. Shifman, N. G. Uraltsev and A. I. Vain- [12] 
shtein, Int. J. Mod. Phys. A9, 2467 (1994); M. Neubert, 
Phys. Rev. D 49, 3392 (1994). [13] 

[4] M. Neubert, Phys. Rev. D49, 4623 (1994). 

[5] A. K. Leibovich, I. Low and I. Z. Rothstein, Phys. Rev. [14] 
D 61, 053006 (2000) [15] 
A. Bornheim et al. [CLEO Collaboration], 



[6] 

[7] G W. Bauer, M7 E. Luke and T. Mannel, 



arXiv : hep-ex/0202019 



arXiv:hep-ph/0 1020891 

R. K. Ellis, W. Furmanski and R. Petronzio, Nucl. Phys 
B 207, 1 (1982), Nucl. Phys. B 212. 29 (1983L 
F. Campanario and T. Mannel, arXiv:hep-ph/0201136 
M. E. Luke and A. V. Manohar, Phys. Lett. B 286, 348 
(1992) 

A. V. Manohar and M. B. Wise, Cambridge Monogr. 

Part. Phys. Nucl. Phys. Cosmol. 10, 1 (2000). 

M. Gremm and A. Kapustin, Phys. Rev. D 55, 6924 

(1997) 

C. Bauer, Phys. Rev. D 57, 5611 (1998); Erratum-ibid. 
D 60, 099907 (1999). 

T. Mannel and M. Neubert, Phys. Rev. D 50, 2037 (1994) 
A. K. Leibovich. Z. Ligeti, M. B. Wise, 



arXiv:hep-ph/0205148 



